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Abstract. Granular media are examined with the focus on polydisperse 
mixtures in the presence of two localized heat-baths. If the two driving 
energies are similar, the large particles prefer to stay in the 'cold' regions 
of the system - as far away from the energy source as possible. If one of the 
temperatures is larger than the other, the cold region is shifted towards the 
colder reservoir; if the temperature of one source is much higher, a strong, 
almost constant temperature gradient builds up between the two reservoirs 
and the large particles are found close to the cold reservoir. Furthermore, 
clustering is observed between the heat reservoirs, if dissipation is strong 
enough. 



1. Introduction 

The segregation of granular materials is an effect of eminent importance for 
industrial operations and has been a subject of research for decades. The 
behavior of powders within the industrial environment, e.g. silos, hoppers, 
conveyor belts or chutes, displays interesting effects - one of them is size 
segregation. Segregation or the mixing properties of granular media are not 
yet completely understood and thus cannot be controlled under all circum- 
stances. For a review of experimental techniques, theoretical approaches, 
and numerical simulations see Refs. [1, 2] and the references therein. 
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A lot of effort has been invested in the understanding of size-segregation 
(see this proceedings for a recent overview of the state of the art). It turns 
out that segregation can be driven by geometric effects, shear, percolation 
and also by a convective motion of the small particles in the system [3] . Seg- 
regation due to convection, in rather dilute, more dynamic systems appears 
to be orders of magnitude faster than segregation due to purely geometrical 
effects in dense, quasi-static situations [4,5]. However, there are still a lot 
of open questions which are subject to current research on model granular 
media [1, 6-9]. 

Most of the segregation phenomena are obtained in the presence of 
gradients in density or temperature. Here we isolate the latter case, i.e. we 
examine the segregation of two species of grains in the presence of local 
heat reservoirs, but in the absence of external forces like e.g. gravity. 

2. The inelastic hard sphere model 

In this study, we use the standard interaction model for the instantaneous 
collisions of particles with radii a,k, and mass m^, with the subscripts k = 
S or L, for small and large particles, respectively. This model accounts 
for dissipation, using the restitution coefficient r, and is introduced and 
discussed in more detail in Refs. [9-13]. The post-collisional velocities v' 
are given in terms of the pre-collisional velocities v by 



v i,2 



(1 + r) m red 
Vl,2 T v n , (1) 

mi,2 



with v n = [(vi — V2) • v n ] n, the normal component of vi — V2, parallel 
to n, the unit vector pointing along the line connecting the centers of the 
colliding particles. The reduced mass is here m re d = mi?ri2/(w-i + %)• If 
two particles collide, their velocities are changed according to Eq. (1). 

If a particle i crosses a line of fixed temperature Tj, its velocity is 
changed in magnitude, but not in direction, according to the rule 

v^ = ±^rV (2) 

with the random thermal velocity ±v j drawn from a Maxwellian velocity 
distribution, and with random sign. (If the '+' sign is used in Eq. (2), a 
large net mass flux occurs when a cluster of particles crosses a line of fixed 
temperature. If the '— ' sign is used the two subsystems have conserved 
particle numbers, but are still coupled via collisions across the boundaries.) 



In 2D one has v i = Jv% + Vy , where v x and v y are the components 
of the thermal velocity vector, each distributed according to a Gaussian 
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distribution. Eq. (2) is also applied to particles after a collision, if the 
particle's center of mass is closer than cll to one of the heat reservoirs j. 

In order to obtain the two random velocities v x and v y , two random 
numbers r\ and r2, homogeneously distributed in the interval [0 : 1] are 



used. With the desired typical thermal velocity v = J2Tj/mi, one has 



V — v 2 lnr2Cos(27rri) , and (3) 



v 



V 



V — v 2 lnr2 sin(2-/rri) , (4) 

using the method by Box and Muller, as described in Ref. [14]. 

If the velocity of the particle would be set to v' = (v x ,v y ), artificial 
peaks in the density at the positions Z\ and Z2 would be observed. This 
artefact was the reason to choose the above described way of thermal cou- 
pling. Note, that our coupling to a reservoir does not guarantee a certain 
temperature in a small volume around the reservoir [15, 16]. It rather ad- 
justs the velocity of all particles that came along the reservoir and touched 
it with their center of mass. Those particles which approach the reservoir 
usually have a lower temperature and thus reduce the mean. As r decreases, 
the reduction increases. Our choice of thermal coupling is arbitrary, how- 
ever, the discussion of thermostating is far from the scope of the present 
paper - therefore, we restrict ourselves to the method introduced here. 

Our method of imposing a temperature gradient does not have any sim- 
ple physical analog, but it does allow us to isolate the effects of the temper- 
ature gradient. More physically feasible energy sources, such as vibrating 
walls perturb the motion more than the method used here. 

3. The event-driven simulation method 

For the simulation of the hard spheres, we use the event-driven algorithm, 
originally introduced by Lubachevsky [17], and applied to the simulation 
of granular media e.g. in Refs. [13]. In these simulations, the particles fol- 
low an undisturbed translational motion until an event occurs. An event 
is either the collision of two particles, the crossing of one particle with the 
boundary of a cell (in the linked-cell structure which is used for algorithmic 
optimization only), or the crossing of one of the lines with fixed temper- 
ature. A particle-particle collision is treated as described in the previous 
section, a cell-boundary crossing has no effect on the particle-motion, and 
the crossing of a fixed temperature location leads to a change of velocity 
according to Eq. (2). Only the particle(s) involved in the last event is (are) 
treated and their next event is computed. In the next step (which does 
not imply a fixed step in time), the next event out of all possible events is 
treated. 
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Figure 1. Snapshots from different simulations with r = 1.0, 0.95, 0.90, and 0.60. The 
heat reservoirs are indicated by horizontal lines; the large particles are displayed in light 
grey, the small in black. 



4. Boundary conditions and system parameters 

The simulations presented here (typical situations are displayed in Fig. 1) 
were performed with iV = 1296 particles in a two-dimensional (2D) box 
of size I = l x = l y = 0.04 m. The box has periodic boundaries, i.e. a 
particle that leaves the volume at the bottom (left), immediately enters 
it at the top (right) and vice-versa. Two particle types (S, L) are used 
with as = a = 2.5 x 10 _4 m and ql = 2a = 5 X 10~ 4 m. In the system 
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with dimensionless size L = l/a = 80, N$ = 1176 particles of the small 
species and Nl = 120 particles of the large species coexist. Note that we 
used a rather arbitrary choice for the calculation of the particle mass (see 
above) so that rriL = 8ms f° r the size ratio used here. The fraction of the 
area which is occupied by the particles, i.e. the total volume fraction, is 
v = ir(N s a s + N L a 2 L )/l 2 « 0.2. 

The hot and cold temperature reservoirs are situated at the vertical 
positions Z2 = ^/a = 15 and Z\ = Z\/a = 55, respectively, and reach 
over the whole width of the system. The temperatures of the reservoirs 
are T\ and T2 — but since no external body force like gravity is involved 
the behavior of the system does not depend on the absolute value of T\ or 
T2; only the ratio T2/T1 is important, and it is varied in the range from 
T2/T1 = 1 to 9 in the following. In Fig. 1 snapshots of simulations with 
T2/T1 = 1 and different values of r are plotted. For small r, clustering is 
observed and the averaging over horizontal slices (data are presented in 
Fig. 2) becomes questionable. 
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Figure 2. Density p (top) and temperature r (bottom), plotted as function of the 
vertical position in the system Z — z/a. Data for small (Left) and large (Right) particles 
are presented for T2/T1 = 1 and different restitution coefficients r as given in the inset. 
The data at z/a = and z/a = 80 are identical due to periodic boundaries. 



6 STEFAN LUDING, OLIVER STRAUSS, AND SEAN MCNAMARA 



ii 



T 2 =T, 
T,=2 T, 

t;=3 T, 
t;=st, 

T,=9 T, 



0.06 
0.05 
0.04 
0.03 
0.02 
0.01 F*«SBgg 







10 20 30 40 50 60 70 

z/a 




20 30 40 50 60 70 

z/a 




1.4 


1 | 


t 2 =t, ; - 


1.2 


T 2 =3T, - . 






1 
0.8 






0.6 


N $ ■■•■ 


%.., i 


0.4 


>^K 


■\.-a \. .J"' ; 


0.2 


ysxxaisft^omKm^^^-g^^^^i^^gi^^ 



10 20 30 40 50 60 70 80 



Figure 3. Density p (top) and temperature r (bottom), plotted as function of the 
vertical position in the system Z — z/a. Data for small (Left) and large (Right) particles 
are presented for r = 0.99 and different temperature ratios T2/T1. 



5. Results 



A quantitative measure of segregation is the partial density of the different 
species. We present pk = n^/A^, the particle number density rik(z) weighed 
by the number of particles N^ of species k in Fig. 2. The strength of seg- 
regation increases with decreasing r, only for small r = 0.60, segregation is 
rather weak due to clustering. Fluctuations in density are correlated to the 
heat reservoirs: Close to a "hot" region in the vicinity of a energy source, 
the density is lower than in the "cold" regions in between. Note that the 
large particles segregate, while the small particles make up a "background 
fluid" with a comparatively small density variation, if dissipation is weak. 
Thus, all particles prefer regions of low temperature, but the large ones are 
attracted by the cold regions more strongly. 

Now, the restitution coefficient is fixed to r = 0.99 and the temperature 
ratio T2/T1 is modified. In Fig. 3 the situation is presented for different 
ratios T2/T1 = 1, 2, 3, 5, and 9. The large particles segregate from the 
"background fluid" made up by the small particles and the quality of seg- 
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regation increases with the magnitude of the temperature gradient. For 
Ti > T\ the large (and heavier) particles can be found close to the colder 
heat reservoir, as different to the situation discussed above. When heavier 
particles have the same temperature as light ones, their mean velocity is 
reduced, so that they can not diffuse away from the cold heat reservoir. 

6. Summary and Conclusion 

In summary, we presented simulations of different size particles in the pres- 
ence of a temperature gradient. If the temperature gradient is due to the 
dissipative nature of the material, the large particles move towards the cold 
regions, as far away from an energy source as possible. If the gradient in 
temperature is externally imposed, most of the large particles prefer to 
move towards the colder heat reservoir. When dissipation is strong enough, 
clustering is observed in the locally driven system. 

Further possible studies involve the prediction of the density and tem- 
perature profiles with kinetic theory - first for the monodisperse and later 
also for the polydisperse case. 
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